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ABSTRACT 

The  stability  of  a  high-speed,  axisymmetric  boundary  layer  is  investigated  using 
secondary  instability  theory  and  direct  numerical  simulation.  Parametric  studies  based 
on  temporal  secondary  instability  theory  identify  subharmonic  secondary  instability  as 
a  likely  path  to  transition  on  a  cylinder  at  Mach  4.5.  The  theoretical  predictions  are 
validated  by  direct  numerical  simulation  of  temporally-evolving  primary  and  secondary 
disturbances  in  an  axisymmetric  boundary-layer  flow.  At  small  amplitudes  of  the 
secondary  disturbance,  predicted  growth  rates  agree  to  several  significant  digits  with 
values  obtained  from  the  spectrally-accurate  solution  of  the  compressible  Navier-Stokes 
equations.  Qualitative  agreement  persists  to  large  amplitudes  of  the  secondary  distur¬ 
bance.  Moderate  transverse  curvature  is  shown  to  significantly  affect  the  growth  rate 
of  axisymmetric  "second  mode"  disturbances,  the  likely  candidates  of  primary  instabil¬ 
ity.  The  influence  of  curvature  on  secondary  instability  is  largely  indirect  but  most 
probably  significant,  through  modulation  of  the  primary  disturbance  amplitude. 
Subharmonic  secondary  instability  is  shown  to  be  predominantly  inviscid  in  nature, 
and  to  account  for  spikes  in  the  Reynolds  stress  components  at  or  near  the  critical 
layer. 
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1.  OVERVIEW 


Renewed  interests  in  the  development  of  high-speed  civil  transport  (HSCT)  aircraft  and  the  hyper¬ 
sonic  National  Aerospace  Plane  (NASP)  have  rekindled  research  in  the  area  of  high-speed  transition 
prediction  and  control.  With  increased  Mach  number  and  higher  llight  altitude,  the  length  of  the  transi¬ 
tional  How  region  increases,  extending  over  large  portions  of  aircraft  surfaces.  In  general,  realistic  drag 
and  surface  healing  estimates  arc  essential  to  the  design  of  high-speed  aircraft  and  depend  on  accurate 
prediction  of  transition  locations  and  extent. 

In  this  paper,  we  combine  secondary  instability  theory  (SIT)  with  direct  numerical  simulation 
(DNS)  to  investigate  the  stability  of  high-speed  boundary-layer  How  along  a  prototypical  "fuselage"; 
i.e.,  a  cylinder.  Specifically,  we  focus  on  the  linear  and  initial  nonlinear  stages  of  forced  transition,  via 
a  subharmonic  secondary  instability  route.  We  find  secondary  instability  theory  and  direct  numerical 
simulation  to  complement  one  another  in  several  ways.  First,  in  the  absence  of  high-speed  experimental 
data,  DNS  provides  a  means  to  substantiate  theory  and  to  establish  the  limits  of  its  validity.  Second, 
within  its  region  of  validity,  secondary  instability  theory  is  preferable  to  computationally-intcnsivc 
numerical  simulation  to  explore  the  parameter  space.  Third,  when  the  nonlinear  stages  of  a  controlled 
instability  process  arc  simulated  numerically,  initial  conditions  derived  from  SIT  provide  a  "jump  start", 
thereby  reducing  significantly  the  integration  times  necessary  to  attain  highly  nonlinear  states. 

In  the  next  section  we  review  some  experimental  and  theoretical  results  which  are  relevant  to  ins¬ 
tability  processes  on  cylinders.  After  briefly  defining  the  geometry,  governing  equations,  and  nondi- 
mensionalization  in  the  third  section,  in  Section  4  we  address,  in  general,  how  transverse  curvature 
influences  instability  mechansims.  As  a  preface  to  the  discussion  of  secondary  instability  theory  in  Sec¬ 
tion  6,  we  discuss  relevant  results  from  temporal  linear  stability  theory  for  compressible  flows  in  Sec¬ 
tion  5.  In  Section  7,  results  of  direct  numerical  simulation  arc  presented  and  arc  compared  with  theory. 
Closing  observations  arc  offered  in  the  final  section. 
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2.  RKVIKW 

Whereas  instability  mechanisms  in  incompressible  boundary  layers  have  been  studied  extensively, 
research  interest  in  the  stability  of  compressible  Hows  has  been  sporadic.  Following  the  pioneering 
work  of  Lees  and  Lin  [1]  in  1946,  there  was  a  very  active  period  in  compressible  stability  research 
which  lasted  into  the  early  1960’s.  After  nearly  two  decades  of  relative  dormancy,  interest  in  this  area 
was  revived  in  the  mid  1980’s  when  national  attention  was  directed  toward  HSCT  and  NASP.  Theoret¬ 
ical  papers  remain  relatively  few  and  experimental  results  are  still  quite  rare.  Although  classical  linear 
stability  theory  has  been  thoroughly  validated  for  incompressible  boundary  layers  (most  notably  |2J),  the 
only  existing  confirmations  of  supersonic  and  hypersonic  stability  theory  ol  which  the  authors  are  aware 
are  those  of  Laufer  and  Vrebalovich  (3|  and  Kendall  [4],  The  former  experiment  investigates  super¬ 
sonic  How  over  a  flat  plate;  the  latter  examines  supersonic  and  hypersonic  How  over  a  flat  plate  and 
hypersonic  llow  over  a  cone.  Although  the  agreement  between  linear  theory  and  these  experiments  is 
imperfect,  it  is  encouraging,  given  the  difficulty  of  making  accurate  measurements  in  a  high-speed  flow 
environment.  The  situation  is  well  summarized  by  Laufer  and  Vrebalovich  13)  w'ho  state:  "In  general, 
the  experiments  confirm  the  basic  assumptions  and  predictions  of  the  existing  stability  theory  and  also 
suggest  the  desirability  of  improvement  in  the  theory'  in  certain  phases  of  the  problem".  With  regard  to 
direct  numerical  simulation,  there  have  been  only  a  handful  of  stability  simulations  for  compressible 
flows  (e.g.  [5)),  and  all  of  these  within  the  past  live  years,  as  reported  by  Kleiscr  and  Zang  in  their 
recent  review  article  16] .  In  general,  there  remains  a  great  need  for  carefully  conducted  experiments  on 
compressible  transition,  well-coordinated  with  theoretical  analyses  and  direct  numerical  simulations. 

Concerning  the  cylinder,  it  is  known  from  the  experiment  of  Browm,  conducted  in  1957  and  only 
recently  published  [71,  that  natural  transition  on  a  cylinder  in  a  subsonic,  axial  flow  occurs  via  the  now 
classical  secondary  instability  route  as  described  by  Herbert  and  coworkers  [8,9],  Brown’s  experiment 
shows  explosive  secondary  instability  emerging  from  axisymmetric,  linite-amplitude  Tollmien- 
Schlichting  (TS)  instability  waves.  The  wavelengths  in  the  streamwisc  and  azimuthal  directions  arc 
roughly  equal  and  the  secondary  instability  is  of  subharmonic  (H)  type  as  evidenced  by  a  staggered 
arrangement  of  A-vorticcs.  More  recently,  Kcgclman  and  Mueller  [10]  have  observed  both  fundamental 
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(K)  type  and  subharmonic  (H)  type  secondary  instabilities  in  forced  transition  experiments  on  an  ogive- 
cylinder  in  a  low-speed  How.  At  low  levels  of  forcing  and  in  natural  transition,  the  staggered  vortex 
arrangement  is  favored,  whereas  at  high  forcing  amplitudes  the  tendency  is  toward  the  aligned  (K) 
arrangement.  Qualitatively,  this  is  in  keeping  with  observations  of  transition  in  flat-plate  boundary 
layers  as  reported  in  [6). 

With  regard  to  (.spatial)  stability  theory  for  incompressible  flo.vs  along  cylinders,  Kao  and  Chow 
[111  show  curvature  to  have  a  stabilizing  trend  on  TS  instability  waves,  characterized  by  a  shift  of  the 
neutral  curve  toward  higher  Reynolds  numbers  in  the  frequency-Reynolds  number  stability  plane. 

Experiments  of  transition  in  supersonic  flows  along  cylinders  are  quite  rare.  Among  these  few 
are  the  ballistic  range  experiments  of  Sciff  and  Short  [12]  and  Whitfield  and  Iannuzzi  [13].  In  the 
former  experiment,  remarkable  shadowgraph  photographs  capture  laminar-turbulent  transition,  and,  in 
one  test  case,  a  turbulent  burst.  Curiously,  in  the  latter  experiment,  transition  data  on  cones  are  com¬ 
pared  to  data  for  a  hollow  cylinder,  which  the  authors  claim  to  be  "equivalent  to  a  flat  plate".  Our 
findings  do  not  support  this  equivalence.  There  have  been,  however,  several  transition  experiments  on 
cones  [4,13-18],  which  arc  particularly  enlightening.  Despite  simplifying  assumptions  which  neglect 
some  effects  of  transverse  curvature,  Mack  [  19 j  obtains  qualitative  agreement  between  linear  theory  and 
an  experiment  by  Stetson  and  coworkers  [14]  which  investigates  Mach  8  How  past  a  sharp,  7-dcgree 
half-angle  cone.  The  frequency  of  the  most  amplified  disturbance  and  the  neutral-curve  obtained 
analytically  by  Mack  match  closely  the  experimentally-observed  values.  However,  observed  and 
predicted  growth  rates  differ  considerably.  Another  experiment  of  particular  interest  is  the  flow  visuali¬ 
zation  experiment  of  Fischer  and  Weinstein  [15]  on  a  shallow  cone  in  a  very  high-speed  flow  (frecs- 
tream  Mach  number  of  20)  in  helium.  Fischer  and  Weinstein  observe  transition  to  originate  near  the 
critical  layer  and  to  propagate  inward  to  the  body  over  a  relatively  long  distance.  For  hypersonic  flow 
over  a  cone,  Kendall  [4]  observes  in  the  boundary  layer  "well-ordered  rope-like  waves  ....  which  persist 
for  relatively  long  distances  and  which  culminate  in  transition".  These  latter  two  results  among  others 
suggest  at  least  two  qualitative  differences  between  low-speed  and  high-speed  transition  mechanisms:  at 
high  speeds  the  transition  region  is  relatively  long;  and  at  high-speeds  transition  originates  near  the  crit- 


ical  layer  rather  than  near  the  wall. 
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With  regard  to  the  theory  of  stability  for  high-speed  Hows  along  cylinders,  Duck  |20j  shows 
transverse  curvature  to  have  a  stabilizing  influence  on  inviscid,  axisymmetric  instability  modes  at  (edge) 
Mach  numbers  of  2.8  and  3.8.  On  the  basis  of  spatial  linear  stability,  Malik  and  Spall  121]  show  a 
similar  tendency  for  axisymmetric  "second-mode"  disturbances  at  Mach  5.0.  However,  they  also  find 
that  "first-mode"  instability  waves  can  be  destabilized  by  increasing  curvature.  These  results  are  corro¬ 
borated  using  temporal  linear  stability  theory  by  Macaraeg  and  Daudpota  ( 22]  and  by  Pruett  and  Strectt 
123], 

Finally,  Ng  and  Erlebachcr  [24]  examine  secondary  instability  mechanisms  in  compressible  flat- 
plate  boundary  layers  for  a  variety  of  edge  Mach  numbers  up  to  4.5  and  find  subharmonic  secondary 
instability  to  be  a  likely  path  to  transition  in  high-speed  flow.  Unlike  the  explosive  nature  of  secondary 
instability  in  incompressible  flow,  secondary  instabilities  in  hypersonic  flows  typically  grow  more  gra¬ 
dually,  with  amplification  rates  only  somewhat  larger  than  that  of  the  most  unstable  primary  distur¬ 
bance. 

These  experimental  and  theoretical  findings  provide  a  tentative  basis  with  regard  to  the  investiga¬ 
tion  of  high-speed  instability  processes  on  cylinders  to  follow.  First,  compressible  linear  stability  theory 
is  applicable,  at  least  qualitatively.  Second,  transverse  curvature  may  have  a  non-ncgligible  influence 
on  stability.  And  third,  instabilities  at  high  speeds  may  differ  fundamentally  from  those  at  low  speeds. 

3.  GEOMETRY,  GOVERNING  EQUATIONS,  AND  NONDIMENSIONALIZATION 
Coordinate  System  and  Governing  Equations 

Consider  the  compressible,  axial  flow  along  a  cylinder  of  radius  /?*,  as  shown  in  Fig.  1.  If 
(i/,v*,w‘)  denote  velocity  components  in  the  cylindrical  coordinate  system  (x*,0,  r*),  and  p*  and  T 
denote  the  fluid  density  and  temperature,  respectively,  then  the  flow  is  governed  by  the  compressible 
Navicr-Stokcs  equations  in  cylindrical  coordinates,  as  given  in  Appendix  E  of  Thompson  [25]. 


Throughout  this  discussion,  dimensional  quantities  arc  denoted  by  superscript  *.  Our  interest  is  in  the 
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stability  of  the  boundary  layer  downstream  of  the  bow  shock.  Accordingly,  we  assume  the  post-shock 
boundary-layer  edge  conditions  to  be  constant.  And,  for  convenience,  we  define  the  wall-normal  coor¬ 
dinate  z*  =  r'  -  R' . 


Nondimensionalization 


To  reduce  the  parameter  space  and  to  gain  physical  insight,  it  is  desirable  to  render  the  governing 
equations  dimensionless.  For  the  problem  at  hand,  it  is  most  convenient  to  scale  flow  quantities  by 
their  respective  edge  values,  denoted  by  subscript  "c",  and  to  scale  lengths  by  a  reference  length  L*  to 
be  defined  shortly.  The  following  parameters  then  appear  in  the  dimensionless  governing  equations: 
Me ,  the  edge  Mach  number;  P*,  the  Prandtl  number;  ReL.,  the  Reynolds  number  based  on  and  y, 
the  ratio  of  specific  heats.  Specifically,  these  are  defined  as  follows: 


CpY 4  „  P  eUeL 

— —  ;  Re,.  =  - ; — 

K,  L  ]Xe 


(1) 


In  eqs.  (1),  C‘  and  C‘  are  respectively  the  specific  heats  at  constant  pressure  and  volume,  and  a'  is  the 
speed  of  sound.  The  ratio  of  specific  heats  and  the  Prandtl  number  are  assumed  here  to  be  constant 
with  values  of  1.4  and  0.7,  respectively.  The  viscosity  (J€*  and  the  thermal  conductivity  k‘  are  each 
assumed  to  vary  with  temperature  according  to  Sutherland's  law;  namely, 


e_=*L=  r^l+c  t)  .  T=r_ 

n;  k-;  t+c  i  ’  re 


(2) 


where  dimensional  temperatures  arc  expressed  in  degrees  Rankine  and  7‘*=198.6<7? . 

There  are  several  possible  choices  of  length  scale  L‘ .  For  consistency  with  the  work  of  others, 
we  use 


where  v*=pe’/p*  is  the  kinematic  viscosity.  An  alternate  length  scale  is  the  boundary-layer 
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displacement  thickness  defined  by 


8*  =]|1 

0  Peue 


8*(1+~r)  =  jrrr\l-^-r-r]dz' 
2R  Jn  R  p,u. 


Equation  (4a)  is  the  standard  definition  for  planar  boundary  layers,  whereas  (4b)  applies  to  thick 
axisymmetric  boundary  layers  [26J.  Note  that  (4b)  involves  the  solution  of  a  quadratic  equation  and 
that,  for  very  large  radius  R\  (4b)  degenerates  to  (4a).  At  Mach  4.5  and  8*/^ *=0.1,  for  example, 
equations  (4a)  and  (4b)  diffci  in  8*  by  about  0.5%. 


4.  EFFECTS  OF  TRANSVERSE  CURVATURE 


Whereas  the  boundary  layer  along  a  fiat  plate  is  self-similar  in  the  absence  of  a  streamwise  pres¬ 
sure  gradient,  that  along  a  cylinder  is  non-similar.  It  can  be  argued  that  self-similarity  is  possible  only 
for  geometries  in  which  there  is  no  inherent  length  scale,  as  is  the  case  for  semi-infinite  fiat  plates  and 
wedges.  In  contrast,  the  cylinder  has  a  characteristic  radius  R ' .  As  the  boundary-layer  displacement 
thickness  8*  grows  in  the  streamwise  direction  while  the  radius  remains  fixed,  their  ratio 


R 


by  which  wc  quantify  the  curvature  C ,  increases  in  the  streamwise  direction,  destroying  self-similarity. 

Transverse  curvature  C  affects  boundary-layer  stability  in  two  ways:  "directly",  through  the  cur¬ 
vature  terms  which  appear  in  the  linearized  disturbance  equations  when  expressed  in  the  natural  cylindr¬ 
ical  coordinate  system;  and  "indirectly",  through  streamwise  evolution  of  the  non-similar  mean  flow. 
Figure  2,  obtained  from  the  spcctrally-accuratc  boundary-layer  code  of  Pruett  and  Sireclt  [23],  shows 
the  streamwise  evolution  of  the  mean  streamwise  velocity  and  the  mean  temperature  in  the  boundary’ 
b<ycr  on  a  cylinder  in  a  high-speed  flow.  The  edge  Mach  number  is  4.5  and  the  edge  temperature  is 
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110°/?  (wind  tunnel  conditions).  Each  curve  of  the  family  corresponds  to  a  unique  curvature  C,  or 
equivalently,  to  a  unique  streamwise  station  x’ .  Thus,  it  is  useful  to  think  of  curvature  as  a  function  of 
the  streamwise  coordinate;  i.e.,  C=C(x‘).  From  cqs.  (3)  and  (5)  we  obtain 


C  = 


V ReR 


i 

.  v  R-  l' 


(6) 


where 


Re 


H  * 


P  ’u’R  ’ 


(7) 


From  eq.  (6)  it  is  dear  that  C(x*)  depends  on  three  dimensionless  quantities;  the  Reynolds  number 
based  on  radius;  the  dimensionless  streamwise  coordinate,  also  based  on  radius  (i.e.,  a*  /?*);  and  die 
ratio  of  length  scales  5*/L*.  Of  these  three  quantities,  the  first  two  are  free  parameters,  whereas,  as  it 
will  become  clear,  the  latter  is  not.  Typically,  for  aerodynamic  applications.  C«l.  However,  small  C 
does  not  imply  immediately  that  curvature  effects  are  negligible,  as  will  be  shown. 

Table  I  helow  compares  the  flight  environments  and  geometric  characteristics  of  a  variety  of 
high-speed  flight  vehicles.  In  the  table,  is  the  cruise  (frecstrcam)  Mach  number,  h  is  the  cruise 
altitude  in  /txlO3,  R'miX  is  the  maximum  radius  of  the  fuselage  in  feet.  Re  x  is  the  frecstrcam  unit  Rey¬ 
nolds  number  at  cruise  conditions  in  units  of  /t~'xl06,  and  A  is  the  fuselage  aspect  ratio  based  on  the 
maximum  radius;  i.e.,  A=t' //?„„,  where  /*  is  the  length  of  the  fuselage. 


Flight  Environments 

Vehicle 

h 

Re  ] 

fima. 

4 

F-15 

1.55 

37 

3.54 

- 

Boeing  HSCT  (27 1 

1.7 

44 

3.01 

'6.7 

'45. 

MD  HSCT  [28] 

3.2 

65-70 

1.50-2.03 

'8.4 

'37. 

NASA  HSCT  (29] 

2.62 

60-68 

1.42-2.17 

5.765 

54.6 

Reentry  Cone  [18] 

20.1-18.2 

140-40 

.3 1  -40. 

1.14 

11.4 

NASP 

>10 

.2-. 5 

_ 

>20 

Table  I.  Flight  environments  of  various  high-speed  vehicles 
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In  general  the  trend  lor  all  three  quantities  whieh  allecl  curvature  is  lowiird  greater  iniluence  with 
increased  Mach  number.  First,  as  the  table  shows,  unit  Reynolds  numbers  tend  to  decrease  with  die 
higher  operational  altitudes  of  taster  vehicles.  Consequently,  Re  .  tends  to  decrease  with  increasing 
Mach  number.  Second,  high-speed  vehicles  tend  to  be  slender,  and  large  aspect  ratio  implies  that  die 
ratio  x'  R'  can  attain  large  values.  Third,  and  perhaps  most  significant,  the  relative  boundary-layer 
displacement  thickness  6’ ,L'  grows  as  A//,  as  shown  in  Fig,  A  for  How  over  a  Hat  plate.  (Figure  A 
further  shows  the  significant  effect  of  on  displacement  thickness  by  comparing  results  based  on  typi¬ 
cal  wind  tunnel  and  Ilight  edge  temperatures.)  For  How  along  a  cylinder,  the  ratio  6"//.*  is  also  a  slow 
function  of  the  curvature  itself  as  shown  in  Fig.  4.  Thus  eq.  <b)  is  a  weakly  nonlinear  equation  lor  C . 
given  Mt. ,  ReK.,  and  .v *  :R  * .  For  a  high-speed  Ilight  vehicle  such  as  the  proposed  NASP,  where  A/  „»l, 
Re „.-()(  10r’).  and  ,4>l,  then  C-OfUT1)  at  aft  fuselage  stations.  In  die  next  section,  primary  instabil- 

A 

ity  growth  rates  will  be  shown  to  be  sensitive  to  curvature  of  this  order. 

5.  PRIMARY  INSTABILITY  THEORY 

In  dus  section,  we  consider  the  stability  of  an  infinitesimally  small  primary  disturbance  superim¬ 
posed  on  a  steady,  laminar,  compressible  mean  flow.  The  mean  velocity  and  temperature  profiles  are 
obtained  from  the  numerical  solution  of  the  2D  or  axisymmetric  compressible  boundary-layer  equations 
by  die  spectrally-accurale  method  of  I2.4|. 

Although  the  mean  (low  is  in  general  a  function  of  both  the  streamwise  and  wall-normal  coordi¬ 
nates,  we  invoke  the  classical  parallel  (low  assumption  whereby  the  streamwise  variation  of  the  mean 
How  is  considered  negligible  relative  to  the  transverse  variation,  a  reasonable  approximation  for 
sufficiently  large  Reynolds  numbers  and  suflicicndy  short  disturbance  wavelengths.  As  a  consequence, 
the  primary  disturbance  can  be  written  in  the  modal  form 

?,((  A.9,2  )  -  q>lU)e,(‘“+t« +  c.c.  (8) 


where 
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</  -  ip,«  .v  ,»v  j\'  (9> 

and  similarly  lor  the  vector  g>.  In  eq.  (S)  u  u*  L’  and  are  stream  wise  and  azimuthal 

wavenumbers,  respectively,  i  i'  /.  *  u’  is  dimensionless  time,  o>  (o*  u'  L "  is  dimensionless  circular  I  re  - 
queue),  $  is  a  complex  function  which  defines  the  structure  ol  the  disturbance  in  the  wall-normal  direc¬ 
tion,  and  "c.c."  denotes  the  complex  conjugate  of  the  preceding  terms.  The  subscript  "1"  reiers  to  "pri¬ 
mary"  instability,  in  contrast  to  "secondary"  instability,  denoted  by  subscript  "2"  in  die  next  section,  l  or 
clarity,  when  no  contusion  is  likely  to  result,  unsubscriptcd  wavenumbers  and  frequency  refer  to  the  pri¬ 
mary  instability  mechanism.  To  preserve  the  azimuthal  periodicity  of  How  along  a  cylinder,  the  product 
(J/y  must  be  an  integer. 

Two  theories  are  commonly  used  to  study  instability  mechanisms.  Spatial  theory  assumes  co  to 
be  real,  while  a  and  (i  are  complex.  The  real  part  of  a  is  the  stream  wise  wave  number,  the  negative  ol 
the  imaginary  part  is  the  spatial  growth  rate.  In  contrast,  in  temporal  theory,  <x  and  P  are  taken  as  real 
while  to  is  complex.  The  real  and  imaginary  parts  of  to  determine,  respectively,  the  temporal  frequency 
and  growth  rate  of  the  wave.  In  temporal  theory,  eq.  (Si  defines  a  wave  uidi  a  real  total  wavenumber 
Ci-  xtrTjP  which  propagates  with  complex  phase  velocity  r  coTx  at  the  angle  tan "'((h a)  to  die 
stream  wise  axis  Temporal  and  spatial  theories  give  equivalent  results  only  for  neutrally -stable  waves. 
An  approximate  relationship  between  temporal  and  spatial  theories  is  given  by  the  Cluster  transformation 
| TO),  which  is  increasingly  valid  as  growth  rates  tend  toward  zero.  In  the  present  analysis,  we  consider 
only  temporally-growing  instabilities. 

The  formulation  of  the  eigenvalue  problem  for  die  analysis  of  temporally-evolving  instability 
waves  in  compressible  Hows  is  well  documented  in  Mack  1  11  |  and  its  numerical  solution  v  ia  spectral 
collocation  methods  is  described  thoroughly  in  Macaraeg  et  al  |  T 2 1 .  All  results  presented  in  this  section 
were  obtained  using  temporal  linear  stability  codes  which  are  variants  ol  dial  described  in  |32|,  to 
which  the  reader  is  referred  for  further  detail. 

A  comprehensive  review  of  the  linear  stability  ol  compressible  boundary  layers  can  be  found  in 
the  w'orks  ol  Mack  (31 1  and  Malik  1 3 3 1  l  or  completeness,  some  relevant  results  are  restated  here.  In 
brief,  there  are  at  least  two  fundamental  ways  in  which  instability  mechanisms  in  incompressible  and 
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compressible  boundary  layers  differ.  First,  unstable  inodes  in  compressible  \s all-bounded  Mow  may 
exhibit  either  viscous  or  inviscid  character.  Second,  lor  high-speed  compressible  How,  multiple  instabil¬ 
ity  modes  may  coexist. 

Neutrally-Stable  Modes 


Unlike  the  incompressible  Blasius  boundary  layer,  whose  velocity  profile  is  noninllectional,  die 
compressible  Mat-plate  boundary  layer  possesses  a  "generalized  inlle'  tion  point"  defined  by 


dz 


)  -  0 


(10) 


where  subscript  "0"  denotes  time-averaged  (mean)  quantities.  As  lirst  observed  by  Lees  and  Lin  |1], 
the  presence  of  a  generalized  inflection  point  is  necessary,  and,  with  additional  mild  restrictions, 
sufficient  to  guarantee  die  existence  of  an  inviscid  neutral  mode.  A  further  generalization  of  condition 
(10)  to  compressible  Mow  along  a  cylinder  can  be  found  in  Duck’s  recent  work  ( 20] .  In  addition,  there 
exist  nomnlleciional  inviscid  neutral  modes.  Inflectional  neutral  modes  travel  with  a  phase  velocity 
equal  to  the  velocity  of  the  mean  Mow  at  the  generalized  inflection  point,  whereas  noninllectional  neu¬ 
tral  modes  travel  with  phase  velocity  equal  to  the  edge  velocity  u'  |34|.  The  practical  implication  for 
linite  ReL .  of  the  existence  of  both  noninllectional  and  inflectional  inviscid  neutral  modes  is  that,  in 
wall-bounded  supersonic  Mow.  unstable  modes  may  exhibit  cither  viscous,  inviscid,  or  mixed  character. 
In  contrast.  Tollmien-Schlichting  (TS)  instability  modes  in  the  incompressible  flat-plate  boundary  layer 
are  purely  viscous  in  origin.  Here  tin  instability  mode  is  considered  to  be  inviscid  in  origin  if  its 
growth  rate  increases  monotonieally  as  Re  .  increases,  whereas  it  is  considered  viscous  if  the  instability 
vanishes  as  Re L. 


The  "First-Mode"  Instability 


In  some  sense,  the  so-called  "first-mode"  instability  is  the  analog  in  compressible  Mow  to  the  TS 


instability  in  incompressible  Mow.  For  edge  Mach  numbers  below  about  1.6,  the  inflection  point  lies 
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close  to  the  wall  and  the  first-mode  instability  is  cf  essentially  viscous  nature.  As  Mach  number 
increases  the  viscous  instability  mechanism  weakens  and  the  first  mode  instability  acquires  invis-  d 
character,  as  the  generalized  inflection  point  migrates  further  from  the  wall.  Above  about  Me=  2.2  invis- 
cid  amplification  rates  increase  dramatically  with  increasing  Me .  For  edge  Mach  numbers  above  about 
3.8,  the  first-mode  instability  can  be  considered  purely  inviscid,  in  the  sense  that  viscosity  has  only  a 
stabilizing  influence  over  the  entire  range  of  a.  Finally,  unlike  TS  instabilities  in  incompressible  flow, 
first-mode  instabilities  in  supersonic  flow  are  most  unstable  when  oblique  (P/0)  |31). 

The  "Second-Mode"  instability 

For  high  supersonic  and  hypersonic  Mach  numbers,  there  exist  high-frequency  modes  of  instabil¬ 
ity  in  addition  to  the  "first"  instability  mode.  The  first  of  the  so-called  "higher"  modes  is  termed  the 
"second"  mode.  In  the  flat  plate  boundary  layer,  it  first  appears  at  about  Me=2.2  and  becomes  more 
unstable  with  increasing  Mach  number  until  about  Me=4.5,  where  it  is  maximally  unstable.  The  second 
mode  dominates  the  first  mode  somewhere  in  an  interval  of  about  A/c=3  to  Afc=4,  depending  upon  edge 
temperature  and  other  conditions.  It  remains  the  most  "dangerous”  (unstable)  mode  for  high  hypersonic 
flow.  Unlike  the  first  mode,  the  second  mode  is  most  unstable  when  it  propagates  in  the  direction  of 
the  mean  flow  r  1  \. 

For  the  flat-plate  boundary  layer  at  Me=4.5,  the  regions  of  instability  for  2D  ((3=0)  first  and 
second  modes  are  distinct,  as  showm  in  Fig.  5,  adapted  from  Mack  |31|.  The  first  and  second-mode 
"neutral -curves"  correspond  to  parameter  pairs  (aJ?eL-)  for  which  disturbances  of  the  form  (8)  are  neu¬ 
trally  stable.  Parameter  values  corresponding  to  points  within  (without)  the  neutral-curves  identify 
exponentially  growing  (decaying)  disturbances.  The  lower  and  upper  sections  of  each  neutral  curve  arc 
commonly  referred  to  as  "branch  1"  and  "branch  2",  respectively.  The  points  on  the  upper  and  lower 
branches  corresponding  to  the  inviscid  limit  (tfe^.— >°°)  identify  respectively  the  "inflectional"  and 
"noninflcciional”  inviscid  neutral  modes  described  in  a  previous  subsection  |34], 
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Eigenfunctions 

Figures  6  compare  and  contrast  the  eigenfunctions  of  highly-amplified  first  and  second-mode  dis¬ 
turbances  in  high-speed  flow  along  a  cylinder.  The  parameters  of  the  flow  are 

Me  =  4.5  ;  ReL.  =  989.4  (Rc&.  =  10000.)  ;  7e*  =  1 10°/?  (11) 

Pr  =0.7  ;  y=1.4  ;  C  =  0.1 

The  first-mode  disturbance  has  A.=58°  and  wavenumbers  (a, p)=(0.0495, 0.0792).  The  second-mode  dis¬ 
turbance  is  axisymmetric  with  wavenumbers  (a,f})=(0.2276,0.0).  it  is  interesting  to  note  that  the  tem¬ 
perature  and  density  components  of  the  eigenfunctions  strongly  dominate  all  other  components,  a  distin¬ 
guishing  characteristic  of  instability  waves  in  high-speed  boundary  layers.  For  the  second  mode  at 
Aft.=4.5,  the  maximum  amplitude  of  the  temperature  fluctuation  is  more  than  15  times  that  of  the 
sircamwise  velocity  fluctuation.  The  amplitude  peaks  of  the  temperature  fluctuations  occur  very  near 
their  respective  critical  layers  (where  phase  velocity  equals  mean  streamwise  velocity),  "\s  denoted  by 
the  vertical  dashed  lines  in  Figs.  6.  This  suggests  that  the  critical  layer  plays  a  special  role  in  high¬ 
speed  transition.  It  is  also  interesting  to  note  that  the  vertical  component  of  the  perturbation  velocity  is 
insignificant  in  the  first-mode  case,  whereas  it  has  magnitude  roughly  equal  to  the  streamwise  perturba¬ 
tion  velocity  in  the  second-mode  case.  The  eigenfunctions  for  the  moderate  curvature  (C=0.1)  case 
shown  in  Figs.  6  differ  only  slightly  from  their  respective  flat-plate  (C=0.0)  counterparts. 

Effects  of  Transverse  Curvature 

We  now  address  how  transverse  curvature  affects  stability,  relative  to  the  stability  characteristics 
of  flow  over  a  flat  plate.  In  Fig.  7,  we  vary  a  to  extract  a  vertical  "slice"  through  Fig.  5  at  a  fixed 
Reynolds  number  .= 955.7,  denoted  in  Fig.  5  by  the  vertical  dashed  line.  The  solid  and  dashed  lines 
of  Fig.  7  contrast  the  results  for  flat-plate  and  cylindrical  geometry  (C=0.1),  respectively.  In  both 
cases,  it  is  clear  that,  at  Mach  4.5,  second-mode  disturbances  arc  the  most  "dangerous",  with  growth 
rates  live  times  their  2D  (axisymmetric)  first-mode  counterparts  and  roughly  twice  that  of  the  most 
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unstable  oblique  (helical)  first  modes.  There  is  a  twofold  effect  of  transverse  curvature  on  axisym- 
metric  first  and  second  modes:  a  stabilizing  trend  characterized  by  a  suppression  of  the  first  and 
second-mode  peaks;  and  a  shift  of  the  most  unstable  wavenumbers  toward  shorter  wavelengths.  In  con¬ 
trast,  the  oblique  first  mode  at  these  conditions  is  destabilized  by  increasing  curvature.  These  trends 
have  been  observed  and  previously  reported  by  Malik  and  Spall  [21),  who  use  spatial  stability  theory. 

In  Figs.  8,  we  examine  the  sensitivities  of  the  growth  rates  of  two  second-mode  disturbances  of 
fixed  wavenumber  to  variations  in  the  curvature  C .  In  principal,  this  is  analogous  to  forcing  the  boun¬ 
dary  layer  at  a  fixed  frequency  and  observing  the  change  in  growth  rate  at  a  fixed  station  as  the  radius 
of  the  cylinder  diminishes.  As  in  Fig.  7,  /?e£. =955.7  and  Me= 4.5.  For  a=0.20,  near  the  lower  branch  of 
Fig.  7  where  the  rate  of  change  of  the  growth  rate  with  respect  to  wavenumber  is  large,  the  variation  in 
C  has  significant  effect  on  the  growth  rate  (Fig.  8a).  In  contrast,  the  effect  of  curvature  is  only 
moderate  for  a=0.22,  near  the  maximally  unstable  wavenumber  (Fig.  8b).  Figures  8  also  show  the  rela¬ 
tive  contributions  of  the  "direct"  and  "indirect"  effects  of  curvature,  as  defined  in  Section  4,  to  the 
"combined"  effect.  To  obtain  the  "direct"  influence  curve,  we  use  a  flat-plate  mean  flow  and  introduce 
curvature  only  in  the  linearized  disturbance  equations.  In  contrast,  the  "indirect"  influence  curve  results 
if  we  account  for  curvature  in  the  mean  flow  but  not  in  the  disturbance  equations.  Qualitatively,  the 
indirect  and  direct  effects  on  the  second  mode  appear  to  be  additive. 

We  close  this  section  with  a  brief  discussion  of  how  temporal  linear  stability  theory  is  adapted  to 
approximately  model  spatially-evolving  disturbances  in  a  growing  boundary  layer.  In  temporal  linear 
theory,  a  disturbance  of  fixed  dimensional  wavelength  (wavenumber)  is  the  approximate  analog  of  a 
fixed-frequency  disturbance  in  spatial  theory.  In  Fig.  5,  each  ray  emanating  from  the  origin  corresponds 
to  a  unique  dimensional  wavenumber,  three  of  which  are  identified  as  a/  <  <  a'h.  Points  along  each 

ray  correspond  to  unique  values  of  the  length  scale  L" ,  or  equivalently  by  eq.  (3),  to  a  unique  stream- 
wise  station  x' .  The  low- wavenumber  (low-frequency)  disturbance,  denoted  by  a',  can  trigger  only  a 
first-mode  instability,  whereas  the  high-wavenumber  disturbance,  denoted  by  a*,  can  give  rise  only  to  a 
second-mode  instability.  However,  the  disturbance  of  moderate  wavenumber  a'm  can  trigger  either  first 
or  second-mode  instabilities,  depending  upon  C . 
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In  a  growing  boundary  layer,  the  increase  in  £'  associated  with  increasing  x’  leads  naturally  to 
the  growth,  stabilization,  and  decay  of  disturbances  of  certain  wavenumbers  a’ .  For  example,  in  Fig.  5, 
point  £  on  the  lower  branch  (branch  1)  identifies  the  value  of x'  beyond  which  a  small -amplitude  dis¬ 
turbance  of  wavenumber  al  becomes  unstable  and  begins  to  grow  exponentially.  At  the  value  of  x * 
corresponding  to  point  £  on  the  upper  branch  (branch  2),  the  disturbance  again  passes  through  neutral 
stability.  Further  downstream  the  disturbance  is  stable  and  its  amplitude  will  decay.  If  nonlinear 
effects  are  insignificant,  the  wave  attains  its  maximum  (presumably  finite)  amplitude  at  point  £,  and  the 
ratio  of  amplitudes  at  points  £  and  £  is  given  by 


Ae 


*f 

j  <0;  (i  )dx 


(12) 


where  w,(x)  is  the  growth  rate  as  a  function  of  position  x  along  the  ray.  Nonlinear  saturation  mechan¬ 
isms  may  dramatically  alter  the  scenario  above  if  wave  amplitude  becomes  large;  however,  because  of 
the  relatively  low  growth  rates  of  instabilities  in  high-speed  boundary -layer  flow,  linear  theory  is  prob¬ 
ably  adequate  to  describe  the  primary  instability  mechanism  unless  the  flow  is  forced  at  very  high 
amplitudes. 

Figs.  5,  7,  and  8,  and  cq.  (12),  taken  together,  imply  that  moderate  transverse  curvature  may 
significantly  alter  the  long  term  temporal  or  spatial  evolution  of  primary  instability  waves,  particularly 
forced  instability  waves.  In  the  next  section  w'C  show  the  onset  of  secondary  instability  to  be  quite  sen¬ 
sitive  to  the  amplitude  of  the  primary  disturbance.  Thus,  curvature  may  exert  an  indirect  yet  significant 
influence  on  transition. 


6.  SECONDARY  INSTABILITY  THEORY  (SIT) 

The  secondary  instability  theory  on  which  the  results  of  this  section  arc  based  is  a  generalization 
to  compressible  flow  of  the  work  of  Herbert  |9|.  Details  can  be  found  in  the  recent  work  of  Ng  and 


F.rlebacher  1 24 1 ,  which  is  summarized  below. 
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We  consider  the  stability  with  respect  to  small  perturbations  of  an  axisymmetric  base  flow 
comprised  of  a  steady,  parallel  mean  flow  onto  which  a  finite-amplitude,  "neutrally-stable"  primary  dis¬ 
turbance  is  superimposed.  Specifically,  following  the  notation  of  the  previous  section,  the  base  flow  can 
be  written 


c?b(tjc,z)  =  5o(z)  +  e,?1(z)ei(cu‘<0l')  +  c.c.  (13) 

where  subscripts  "0"  and  ”b"  refer  to  ''mean'1  and  "base”  flows  respectively,  and  "e/'  is  the  amplitude  of 
the  primary  instability  wave  (normalized  so  that  the  maximum  amplitude  of  the  temperature  fluctuation 
is  unity).  Here  we  consider  only  axisymmetric  primary  waves,  although  in  the  more  general  theory 
[24],  the  primary  wave  can  be  oblique  (helical).  Again  we  restrict  consideration  to  the  temporal  model 
of  stability  in  which  the  wavenumber  a  is  real.  In  a  frame  of  reference  x'  -  x  -  Cjt  which  moves  with 
the  phase  velocity  c  ,=Re(a)i/a)  of  the  primary  wave,  the  base  flow  is  independent  of  lime  and  periodic 
in  the  streamwise  spatial  dimension.  Consequendy,  the  stability  analysis  is  amenable  to  Floquet  theory 
[9],  which  applies  to  differential  equations  with  periodic  coefficients.  With  guidance  from  Floquet 
theory,  a  general  ansatz  for  a  secondary  disturbance  ensemble  of  N  wave  components  is  as  follows: 


^2(f r* ,6,z )  =  e  ,0i2' e  2  e,|W8  £  <j>2 ,>(z)e'/c“  +  c.c.  ( N  even ) 

j=~N/2 


(14a) 


^2(/yC,0,z)  =  e  LU>2‘ e  2  e 


ipse 


(.N-m 

)=-(AM)/2 


+  C.C. 


(. N  odd ) 


(14b) 


where  eqs.  (14a)  and  (14b)  are  valid  for  an  even  and  odd  number  of  wave  components  respectively.  In 
eqs.  (14),  subscript  "2"  denotes  "secondary"  disturbance,  and  h  is  the  "detuning  parameter".  Hen¬ 
ceforth,  p  without  subscripts  refers  to  the  azimuthal  wavenumber  of  the  secondary  disturbance.  For 
h=l,  the  secondary  disturbance  is  of  subharmonic  (H)  type,  whereas  for  h=  0  it  is  of  fundamental  (K) 
type.  For  0  <  h  <1,  eqs.  (14)  model  so-called  "detuned"  secondary  instability  modes.  In  practice,  the 
parameters  of  the  neutrally-stable  primary  wave  correspond  to  some  point  along  the  upper  branch  of  the 
neutral-stability  curve,  in  keeping  with  the  discussion  of  Fig.  5  in  the  previous  section. 


Let  us  now  choose  wavenumbers  a  and  (3  and  suppose  If  -  if b  &  z2Ifi,  where  tfb  and  If2  arc 
defined  in  eqs.  (13)  and  (14),  respectively,  and  where  e2  *s  the  amplitude  of  the  secondary'  disturbance, 
presumed  small.  Substitution  of  If  into  the  compressible  Navier-Stokes  equations,  neglect  of  nonlinear 
terms  in  ef  and  e2,  and  discretization  in  the  z  direction,  results  in  a  complex  generalized  eigenvalue 
problem  tor  the  eigenvalue  co2  and  the  discrete  eigenvector  (j>2  =  l<i>2,o>  <j>2.-i»  ?2,o  l7  >  namely 

(A  3 D  ^  +  A  2D  ^  +  A  i)<J)2  —  tovl  y<J)2  (15) 

In  eq.  (15),  Dp  denotes  a  p  -laid  discrete  spatial  differentiation  operator,  and  the  A,  ’s  are  complex 
malices  of  dimension  5MN  x  5 MN ,  where  M  is  the  number  of  discrete  grid  points  in  z .  The  nearly 
one  thousand  lines  of  FORTRAN  code  required  to  evaluate  the  coefficients  of  the  A ,  matrices  arc  gen¬ 
erated  with  the  help  of  MACSYMA,  a  symbolic  manipulation  language,  without  which  the  analysis  of 
compressible  secondary  instability  would  be  a  formidable  task.  Even  with  a  highly-accurate  numerical 
method,  such  as  the  spectral  collocation  method  in  current  use,  in  which  typically  Af  =73,  the  memory 
requirement  for  an  eigenvalue  calculation  of  N  >6  becomes  prohibitively  large.  Fortunately,  only  a  very 
few  components  of  (14)  are  of  major  importance,  as  suggested  in  Table  II. 


Eigenvalue  vs.  Number  of  Components 

N 

Im(to2) 

-Rc(to2) 

2 

0.002416 

<ur5 

3 

0.002411 

4 

0.002406 

<  1  o_:> 

Table  II.  Variation  of  the  Eigenvalue  of  a  Subharmonic  Secondary  Instability 
with  Number  of  Components  N . 


The  results  presented  in  Table  II  are  for  subharmonic  (h=l)  secondary  instability.  The  flow  parameters 
arc  the  same  as  in  eqs.  (11).  The  primary  disturbance  is  an  axisymmetric  second  mode  of  wavenumber 
a=0.2523,  amplitude  £[=.085,  and  frequency  C0j  =  0.2275  +  / 0.436x10^.  The  azimuthal  wavenumber 
of  the  secondary  disturbance  is  P=0.2673.  Due  to  memory  limitations  and  contamination  by  round-off 


error,  it  is  not  computationally  feasible  to  include  enough  components  in  the  eigenvalue  problem  to 
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demonstrate  convergence  of  the  ansatz  (14)  in  a  rigorous  mathematical  sense.  Rather,  we  have  verified 
the  eigenvalues  by  comparison  against  spectrally-accurate  direct  numerical  simulation.  With  the  inclu¬ 
sion  of  just  two  (N-2)  wave  components  in  the  secondary  disturbance,  the  eigenvalue  is  correct  to  two 
significant  digits.  That  Re(co2)  is  essentially  zero  in  the  frame  x'  implies  that  the  primary  and  secon¬ 
dary  disturbances  are  phase  locked. 

Figure  9  displays  a  rather  typical  scenario  for  secondary  instabilities  at  subsonic  and  low  super¬ 
sonic  Mach  numbers.  The  parameter  values  for  this  test  case,  in  which  four  components  (N- 4)  have 
been  retained  in  the  ansatz  (14),  are 

Me  =  1.6  ;  ReL-  =  613.2  (Reg.  =  1641)  ;  T’e  =  520°/?  (16) 

oc  =  0.0867  ;  p  =  0.1336  ;  C  =  0.051 

The  subharmonic  secondary  instability  mechanism  is  favored  at  low  to  moderate  amplitudes  of  the  pri¬ 
mary  disturbance,  whereas  the  fundamental  secondary  instability  is  more  unstable  at  higher  primary 
amplitudes.  This  suggests  that  subharmonic  instability  is  more  "dangerous"  in  natural  transition  in  this 
speed  range,  in  keeping  with  experimental  observations  for  low-speed  How.  At  higher  Mach  numbers, 
there  is  some  evidence  that  subharmonic  secondary  instability  is  yet  more  strongly  favored.  Figure  10 
displays  the  growth  rate  of  the  secondary  instability  as  a  function  of  the  detuning  parameter  h  for  a 
fixed,  moderate-amplitude  (e^O.085)  primary  disturbance.  The  parameter  values  of  the  flow  arc  the 
same  as  those  in  eqs.  (11)  except  that  ReL.  =994.6.  The  disturbance  parameters  are  a=0.2506,  PR  =21, 
and  N=2.  Of  all  possible  secondary'  instabilities,  the  subharmonic  type  (h= 1)  is  most  unstable.  More¬ 
over,  no  fundamental  type  instability  exists  for  these  parameter  values.  Consequently,  we  consider  hen¬ 
ceforth  only  secondary  instabilities  of  subharmonic  type. 

The  nature  of  secondary  instability  is  that  it  is  triggered  only  after  the  primary  has  achieved  some 
threshold  amplitude.  This  is  clearly  shown  in  Fig.  1 1  which  presents  the  growth  rate  of  a  subharmonic 
secondary  disturbance  of  N  =2  vs.  primary  wave  amplitude  E).  Again,  the  flow  parameters  arc  the  same 
as  given  in  eqs.  (11),  except  that  ReL*=955.7.  In  this  case,  the  disturbance  parameters  arc  a=0.2400 
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and  (3=0.2000.  For  ex  <  2.5%,  no  secondary  instability  is  present.  For  el  >  2.5%,  secondary  instabili¬ 
ties  exist  which,  at  large  eb  attain  growth  rates  slightly  higher  than  that  of  the  maximally  unstable 
second-mode  disturbance  at  the  same  ReL . .  Figure  11  also  compares  the  results  for  the  flat  plate  with 
those  of  a  cylinder  of  curvature  C=0.1.  It  appears  that  secondary  instabilities  per  se  are  little  affected 
by  moderate  transverse  curvature.  This  is  further  substantiated  in  Fig.  12,  which  shows  only  slight 
influence  of  transverse  curvature  on  the  growth  rate  of  the  secondary'  instability,  for  a  fixed  8.5%  pri¬ 
mary  disturbance  amplitude.  However,  since  the  growth  rate  of  the  secondary  disturbance  is  strongly 
dependent  on  primary  disturbance  amplitude,  which  is  itself  sensitive  to  curvature  effects  (refer  to  Fig. 
5,  eq.  (12),  and  Fig.  11),  the  indirect  influence  of  curvature  on  secondary  instability  may  be  substantial. 
Except  for  the  variation  in  C  and  the  fixed  £|,  the  parameter  values  for  Fig.  12  are  the  same  as  those 
for  Fig.  1 1 . 

It  is  interesting  to  compare  the  eigenfunctions  of  primary  (Figs.  6)  and  secondary  disturbances, 
the  latter  shown  in  Fig.  13.  Here,  the  parameter  values  are  the  same  as  those  given  for  Table  II,  with 
N- 2  and  £j  =0.085.  Like  the  eigenfunction  of  the  primary  disturbance,  the  eigenfunction  of  the  secon¬ 
dary  disturbance  is  dominated  by  the  temperature  component  which  peaks  in  amplitude  near  the  critical 
layer  (z*  =  1.088*).  However,  unlike  the  primary  mode,  the  secondary  eigenfunction  has  very  little 
structure  near  the  wall.  Remarkably,  the  structure  of  the  eigenfunction  of  the  subharmonic  secondary 
disturbance  is  qualitatively  very  similar  to  that  of  a  helical  first-mode  disturbance  (Fig.  6a). 

The  amplitude  of  the  secondary  disturbance  eigenfunction,  peaked  well  away  from  the  wall,  its 
relative  insensitivity  to  changes  in  the  mean  flow  (refer  to  Fig.  12),  and  its  relatively  large  growth  rate, 
suggest  that  the  secondary  instability  mechanism  is  predominantly  inviscid.  This  suspicion  is  confirmed 
by  Fig.  14.  For  a  fixed  primary  disturbance  with  ReL*=955.1,  a=0.2400,  and  £i=0.085,  the  Reynolds 
number  of  the  cigensystem  for  a  secondary  disturbance  of  (3=0.2000  and  N=  2  is  allowed  to  vary 
artificially  from  about  103  to  1.5xl06,  while  all  other  parameters  remain  fixed.  The  growth  rate  of  the 
secondary  mode  increases  monotonically  with  increasing  ReL-  and  asymptotes  at  iaige  tieL.  to  approxi¬ 
mate  values  of  0.0053  for  C=0.0  and  0.0052  for  (7=0.1.  The  figure  confirms  once  again,  that  for  fixed 
primary  disturbance  amplitude,  the  secondary  instability  mechanism  is  minimally  influenced  by 
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moderate  transverse  curvature. 

7.  DIRECT  NUMERICAL  SIMULATION  (DNS) 

The  3D,  compressible  Navier-Stokes  equations  in  primitive  variables  |p,u,v,w,p]  (where  p  is  the 
pressure)  are  solved  by  a  highly-accurate  spectral  collocation  method  developed  to  numerically  study 
instabilities  in  the  flat-plate  boundary  layer  [5].  The  algorithm,  adapted  to  the  cylindrical  (x,0,z)- 
coordinate  system  of  Fig.  1,  treats  either  flat-plate  or  cylindrical  geometry,  depending  on  the  value  of  a 
switch  parameter.  A  temporal  model  of  stability  is  assumed,  in  keeping  with  the  linear  analyses  of  pri¬ 
mary  and  secondary  instabilities  presented  earlier.  As  before,  the  mean  flow  is  regarded  as  locally 
parallel  and  subject  to  temporally-evolving,  spatially-periodic  perturbations.  Since  the  parallel  mean 
flow  does  not  exactly  satisfy  the  viscous  Navier-Stokes  equations,  small  forcing  terms  are  added  to  the 
governing  equations  sufficient  to  suppress  temporal  evolution  of  the  mean  flow.  The  spectral  colloca¬ 
tion  method  employs  Fourier  basis  functions  in  the  periodic  sireamwise-*  and  azimuthal-0  directions, 
and  Chebyshev  polynomial  basis  functions  in  the  aperiodic  wall-normaf-z  direction.  Stretching  and 
clustering  mappings  in  the  (finite)  wall-normal  direction  concentrate  collocation  points  in  regions  of 
severe  gradients;  i.e.,  at  the  wall  and  at  the  critical  layer.  These  mappings  are  detailed  in  [5 J. 

The  governing  equations  and  boundary  conditions  define  an  initial  value  problem  which  is 
integrated  in  time  explicitly,  by  means  of  a  3rd-order,  low-storage,  Runge-Kutta  scheme  [35].  Initial 
conditions  arc  formed  by  superposition  of  a  parallel  mean  flow,  derived  from  a  spcclrally-accurate 
boundary-layer  code  [231,  and  eigenfunctions  of  primary  and  secondary  disturbances,  generated  from 
the  temporal  stability  code  of  [24|,  also  spectrally  accurate.  The  mean  boundary-layer  profiles  are  gen¬ 
erated  assuming  an  adiabatic  wall.  In  the  direct  simulation,  the  temperature  at  the  wall  is  kept  fixed  at 
its  mean  value,  and  no-slip  conditions  arc  imposed  on  the  velocity  components.  Either  Dirichlet  or 
Neumann  conditions  are  enforced  at  the  far-field  boundary,  with  little  effect  on  the  solution  provided 
the  boundary  is  at  sufficiently  large  z.  The  values  of  density  at  the  wall  and  at  the  far-field  boundary 
arc  obtained  by  projecting  forward  in  time  using  the  Runge-Kutta  integrator.  Pressure  at  the  wall  is 
extracted  from  the  equation  of  state. 
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Test  Case 

The  parameters  of  the  test  case  arc  as  follows: 

Me  =  4.5  ;  ReL.=  994.6  ;  T’e  =  1 10°/?  (17) 

Pr  =0.7  ;  y=  1.4  ;  C  =  0.1 

The  Mach  4.5  case  at  a  typical  wind-tunnel  edge  temperature  was  selected  for  several  reasons.  As 
mentioned  previously,  the  second  mode  is  most  unstable  at  Mach  4.5  so  that  the  most  likely  type  of  pri¬ 
mary  instability  is  well-defined.  Furthermore,  this  case  has  been  investigated  previously  for  flow  over  a 

flat  plate  by  Mack  (31],  using  linear  stability  theory,  and  by  Erlebacher  and  Hussaini  [5],  using  direct 

numerical  simulation.  Thus,  there  is  data  against  which  to  compare  the  results  for  cylindrical  flow. 
Finally,  we  assume  a  moderate  degree  of  curvature  (C=0.1),  large  enough  to  influence  stability,  yet 
within  the  range  of  values  encountered  in  high-speed  flight. 

Initial  Conditions 

The  wavenumbers,  amplitudes,  and  eigenvalues  of  the  instability  waves  imposed  initially  arc 
given  in  Table  III  below. 


Initial  Conditions 

Instability 

a 

P 

m 

e 

Primary  (1,0) 

0.2506 

0.0 

0.2261  +  ().00032i 

0.0850 

Secondary  (1/2+1) 

0.1253 

±0.2079 

0.2261  +  0.00253i 

0.0085 

Secondary  (3/2+1) 

0.3759 

±0.2079 

0.2261  +  0.00253i 

0.0085 

Table  III.  Disturbance  initial  conditions 

For  simplicity,  we  adopt  a  modal  notation  whereby  "(/,/)"  refers  to  a  wave  of  wavenumbers  (taj(3). 
The  primary  (1,0)  wave  is  an  axisymmetric  second  mode  of  nearly  neutral  stability,  whose  parameters 
correspond  to  a  point  near  the  upper  branch  of  the  second-mode  neutral  curve.  The  secondary  dislur- 
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bance,  obtained  in  the  manner  of  the  previous  section,  is  of  subharmonic  type  with  p/?  =2 1  and  N= 4. 
The  (1/2,1)  component  of  the  secondary'  wave  ensemble  has  ^60°.  For  this  value  of  X,  the  secondary 
disturbance  is  very  nearly  maximally  unstable  with  respect  to  p.  The  initial  amplitude  e2  of  the  secon¬ 
dary  wave  ensemble  is  ten  times  smaller  than  that  of  the  primary  wave,  based  on  the  relative  peak 
amplitudes  of  the  temperature  fluctuations  of  the  (1,0)  and  (1/2,1)  waves. 

Results 

The  results  shown  here  are  for  a  simulation  of  12x12x65  resolution  with  equally-spaced  points  in 
the  periodic  dimensions  and  with  the  near-wall  collocation  points  in  the  wall-normal  direction  distri¬ 
buted  as  shown  by  the  symbols  in  Fig.  16.  The  bounds  of  the  computational  domain  are 
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* 
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The  fully-explicit  time  advancement  results  in  a  severe  constraint  on  the  time  step  when  used  in  con¬ 
junction  with  the  Chebyshev  spectral  collocation  method.  About  8000  time  steps  per  period  of  oscilla¬ 
tion  of  the  primary  are  required  to  maintain  numerical  stability.  The  simulation  requires  about  38 
minutes  of  CPU  time  on  a  Cray  Y-MP  per  period  of  oscillation.  Resolution  (decay  of  the  spectral 
coefficients)  was  monitored  throughout  the  compulation,  which  was  halted  at  about  38  periods  (24  CPU 
hours)  when  resolution  began  to  deteriorate  rapidly.  Work  is  in  progress  to  continue  the  simulation 
further  into  the  nonlinear  regime  using  a  refined  mesh  and  a  high-order  compact-differencing  scheme  in 
the  wall-normal  direction. 

Fig.  15  shows  the  temporal  evolution  of  the  "energy"  in  various  harmonics  of  the  primary  and 
secondary  disturbances,  based  on  a  Fourier  decomposition  in  the  periodic  dimensions  and  an  integrated 
average  in  the  wall-normal  dimension;  namely. 
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,  2  max 

E„  =  T  V  J  Pulley  I2  +  K,  I2  +  \W„  l2| ~dz  (19) 

where  upper  case  UtJ  represents  the  (/  ,j  )-th  complex  Fourier  component  of  the  streamwise  velocity,  etc. 
Solid  and  dashed  lines  show  results  of  the  simulation  whereas  circles  denote  the  predicted  values  based 
on  linear  theory  for  primary  and  secondary  instabilities.  The  growth  rate  of  the  primary  (1,0)  wave  is 
in  close  agreement  with  theory  for  more  than  30  periods  of  oscillation.  Initially,  the  (1/2,1)  and  (1/2, -1) 
(coincident)  components  of  the  secondary  disturbance  grow-  at  a  rate  w-hich  agrees  to  several  significant 
digits;  with  the  values  predicted  by  secondary  instability  theory.  These  modes  begin  to  depart  toward 
slightly  higher  growth  r  ites  at  around  10-15  periods  of  oscillation  of  the  primary.  Over  the  lengthy 
integration,  the  slightly  non-neutral  primary  increases  in  amplitude  by  a  factor  of  about  1.5  as  shown  in 
F  ig.  16.  The  squares  in  Fig.  15  denote  adjustments  to  the  theory  to  incorporate  the  gradual  increase  in 
amplitude  of  the  primary,  and  arc  obtained  by  integrating  the  following  set  of  ordinary  differential 
equations  for  the  primary  and  secondary  amplitudes  /t  j  and  A2 ,  respectively: 

ci4, 

— —  =  Im(c0iVt,  (20) 

at 

dA  2 

— —  =  lm({02(^i))^2 
at 

In  eqs.  (20),  lm(M|)  is  assumed  constant,  but  Im((02(/1  [))  is  interpolated  from  the  dam  in  Fig.  11.  This 
correction  slightly  overpredicts  the  upw-ard  deflection  in  the  growth  of  the  (1/2,1)  mode.  In  contrast,  the 
(3/2,1)  and  (coincident)  (3/2,- 1)  components  of  the  secondary  disturbance  grow  consistently  at  the  rate 
Im(ci)|  +  C0i),  denoted  by  triangles  in  Fig.  15.  This  is  precisely  the  rate  one  would  anticipate  based  on  a 
weak  quadratic  nonlinear  interaction  between  the  primary  (1,0)  wave  and  secondary  (1/2,1)  wave. 

Qualitative  agreement  between  theory  and  simulation  persists  until  about  period  35,  when  the 
growth  rate  of  the  primary  disturbance  deviates  significantly  from  the  value  predicted  by  linear  theory. 
At  this  stage,  modes  generated  by  nonlinear  interactions  are  attaining  significant  energies,  and  transition 
appears  to  be  imminent.  In  particular,  the  (1,2)  mode,  which  originates  from  the  self-interaction  of  the 
(1/2,1)  wave,  has  an  energy  content  nearly  equal  to  that  of  the  primary  instability  wave.  By  the  time 


the  computation  is  terminated  at  period  37.7,  the  temperature  component  of  the  secondary  (1/2,1)  wave 
has  grown  to  more  than  15  times  its  initial  amplitude,  as  shown  in  Fig.  17. 


Reynolds  Stresses 

Figure  18,  obtained  from  the  results  of  DNS,  presents  the  four  significant  components  of  the  six 
possible  Reynolds  stress  components.  These  are  computed  in  temporal  theory  according  to 
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where  primes  denote  perturbation  variables;  eg.,  u'  =  u-u0.  The  individual  curves  in  the  four  plots 
represent  the  wall-normal  distribution  of  Reynolds  stress  at  intervals  of  6  periods  in  time,  beginning 
with  the  initial  condition  and  ending  at  period  36.  With  the  exception  of  the  Tl  t  component,  all  Rey¬ 
nolds  stress  components  are  peaked  at  the  critical  layer,  denoted  by  dashed  lines.  In  Fig.  19,  we  exam¬ 
ine  die  relative  contributions  of  the  primary  and  secondary  disturbances  to  the  xtti  component.  It  is 
clear  that  the  large  growth  in  Reynolds  stress  at  the  critical  layer  is  due  almost  entirely  to  the  secondary' 
instability  mechanism.  Fig.  19  further  compares  the  results  of  SIT  and  DNS  w  ith  regard  to  Tz  ; . 
Theory  and  simulation  can  be  considered  in  qualitative  agreement  to  about  period  18.  Beyond  that, 
there  is  gradual  broadening  of  the  stress  peak  in  the  DNS  results,  due  presumably  to  nonlinear  effects. 


In  Fig.  20  we  make  use  of  the  Caster  transformation  |3()|,  based  on  the  phase  speed  ('i=Rc(o)i)/a 
of  the  primary  wave,  to  convert  the  temporal  evolution  of  T,  z  to  an  approximate  spatial  distribution. 
For  clarity  in  Fig.  20,  intervals  between  tick  marks  along  the  streamwise  x  axis  arc  10  displacement 
thicknesses  in  length  whereas  intervals  in  the  radial  direction  each  denote  a  single  displacement  thick¬ 
ness.  Moreover,  it  should  be  noted  that  the  contour  intervals  are  different  for  Reynolds  stress  levels 
below  and  above  O.fMxlO"4.  It  is  interesting  to  view  Figs.  19  and  20  in  light  of  the  high-speed  experi¬ 
ment  of  Fischer  and  Weinstein  |!5|,  in  which  transition  is  observed  to  originate  near  the  criucai  layer 
and  to  propagate  inwardly  to  the  body.  These  figures  suggest  a  similar  evolution  to  transition  and  a 
possible  link  between  secondary  instability  theory  and  observed  phenomena. 
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The  initial  and  final  times  of  the  temporal  simulation  correspond  via  the  Caster  transformation 
[30)  to  approximate  initial  and  final  sireamwise  stations  of  x’o  =  \0R'  and  Jtj-19/?'  respectively.  Over 
this  spatial  extent,  the  boundary  layer  would  have  thickened  by  about  35%,  a  change  of  length  scale 
large  enough  to  significantly  alter  the  growth  of  the  primary  disturbance  according  to  Fig.  5.  It  appears 
that  modeling  the  growth  of  the  boundary- layer  is  a  necessity  if  direct  simulation  of  high-speed  transi¬ 
tion  is  to  be  quantitatively  valid.  Within  the  context  of  temporal  DNS,  this  can  be  approximated 
through  time  evolution  of  the  mean  (low.  Alternately,  one  must  resort  to  more  computationally- 
intensive  spatial  DNS. 

Fluid  Rotation 

It  has  been  suggested  by  Morkovin  [34),  that  the  "angular  momenta"  vector  poi  is  more  physically 
meaningful  for  compressible  flows  than  is  the  vorticity  vector  GJ=VxiF.  This  can  be  inferred  as  well 
from  the  presence  of  density  in  the  generalized  inflection  point  condition  (10).  In  Figs.  21  we  compare 

L*m  (iu  (iw 

contours  of  azimuthal  vorticity  toe  =  -^(-r - r— )  with  contours  of  azimuthal  angular  momentum 

6  oz  dx 

pa)e  at  period  37.7  of  the  numerical  simulation.  The  contour  patterns  arc  quite  different,  due  in  part  to 
the  decreased  mean  density  at  the  wall.  Particularly  in  Fig.  21b,  it  is  clear  that  a  highly  rotational  layer 
has  formed  at  the  critical  layer,  with  rotational  peaks  at  intervals  of  the  subharmonic  wavelength  Anla. 

8.  CONCLUSIONS 

Because  of  the  computational  expense  of  direct  numerical  simulation,  our  work  to  date  has  been 
restricted  to  very  few  parameter  sets  and  is  very  much  "in  progress".  Accordingly,  these  comments 
should  perhaps  be  interpreted  as  "observations"  rather  than  as  "conclusions". 

1)  Secondary  instability  of  subharmonic  type  exists  in  high-speed  boundary-layer  flows  over  flat 
plates  and  cylinders,  and  is  a  likely  path  to  transition.  However,  due  to  a  general  stabilizing  trend 
with  increasing  Me ,  secondary  instabilities  at  high-speeds  arc  not  as  "explosive"  as  their  counter- 
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parts  in  incompressible  flow.  The  regions  of  linear  and  nonlinear  growth  are  relatively  long,  and 
the  boundary  layer  may  experience  significant  growth  during  laminar-turbulent  transition.  Growth 
of  the  boundary  layer  needs  to  be  modeled  if  the  results  of  direct  numerical  simulation  are  to  be 
quantitatively  valid. 

2)  For  primary  and  secondary  instability  mechanisms,  the  temperature  and  density  components  of  the 
eigenfunctions  dominate  at  high-speeds.  For  example,  at  Mach  4.5,  axisymmctric  second-mode 
disturbances  exhibit  peak  temperature-fluctuation  amplitudes  which  are  15  times  that  of  stream- 
wise  velocity  fluctuations. 

3)  For  primary  and  secondary  instabilities,  the  peak  amplitudes  of  the  temperature  fluctuations  occur 
at  or  very  near  the  critical  layer.  Thus,  the  critical  layer  plays  an  important  role  in  the  instability 
mechanisms  which  arise  in  high-speed,  wall-bounded  flows. 

4)  The  principal  effect  of  transverse  curvature  on  primary  and  secondary  instability  mechanisms 
seems  to  be  quantitative  rather  than  qualitative.  In  general,  axisymmctric  primary  instability 
modes  arc  stabilized  by  moderate  curvature,  which  in  turn  should  delay  the  onset  of  secondary 
instability.  Curvature  effects  tend  to  become  more  significant  with  decreasing  Reynolds  number 
and  increasing  Mach  number. 

5)  The  amplitude  of  the  eigenfunction  of  a  secondary  disturbance,  with  little  structure  at  the  wall; 

the  relative  insensitivity  of  secondary  disturbances  to  changes  in  the  mean  How;  the  relatively 
large  growth  rates  of  secondary  instabilities;  and  the  persistence  of  the  secondary  instability 
mechanism  as  all  suggest  that  the  secondary  instability  mechanism  is  predominantly 

inviscid. 

6)  In  high-speed  boundary-layer  flow,  secondary  instability  is  responsible  for  large  Reynolds  stresses 
near  the  critical  layer,  and  offers  a  possible  explanation  for  the  phenomenon  observed  experimen¬ 
tally  by  Fischer  and  Weinstein  115].  Moreover,  insight  gained  from  secondary  instability  theory 
may  provide  a  means  to  refine  transition  models  based  on  Reynolds  stresses. 
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7)  Secondary  instability  theory  offers  a  possible  means  of  refining  the  eN  method  for  transition  pred¬ 
iction,  based  on  a  primary  disturbance  amplitude  criterion  rather  than  on  an  amplitude-ratio  cri¬ 
terion. 

We  conclude  by  reiterating  that  there  is  great  need  for  carefully-conducted  compressible  stability 
experiments,  well-coordinated  with  analyses  and  direct  numerical  simulation. 


ACKNOWLEDGMENTS 

The  authors  wish  to  acknowledge  the  invaluable  contributions  of  Drs.  Tom  Zang  and  Craig 
Streett  to  this  effort.  We  also  wish  to  thank  Dr.  Michele  Macaraeg  for  providing  the  compressible 
Navier-Stokes  equations  in  a  form  suitable  for  symbolic  manipulation  using  MACYSMA.  The  first 
author  is  grateful  to  the  National  Research  Council  for  the  Associateship  which  has  made  this  work  pos¬ 
sible. 


-27- 


REFERENCES 


1)  Lees,  L.,  and  C.  C.  Lin,  "Investigation  of  the  Stability  of  the  Laminar  Boundary  Layer  in 
Compressible  Fluid",  NACA  TN-11  IS,  1946. 

2)  Schubauer,  G.  B.,  and  H.  K.  Skrainstad,  "Laminar  Boundary  Layer  Oscillations  and  Stability  of 
Laminar  Flow",  J.  Aeronaut.  Sci.,  Vol.  14,  pp.  69-78,  1947. 

3)  Laufer,  J.,  and  T.  Vrebalovich,  "Stability  and  Transition  of  a  Supersonic  Laminar  Boundary  Layer 
on  an  Insulated  Flat  Plate",  J.  Fluid  Mech.,  Vol.  9,  pp.  257-299,  1960. 

4)  Kendall,  J.  M.,  Jr.,  "Wind  Tunnel  Experiments  Relating  to  Supersonic  and  Hypersonic  Boundary- 
Layer  Transition",  A1AA  J .,  Vol.  13,  No.  3,  pp.  290-299,  1975. 

5)  Eriebacher,  G.,  and  M.  Y.  Hussaini,  "Numerical  Experiments  in  Supersonic  Boundary-Layer  Sta¬ 
bility",  Phys.  Fluids  A,  Vol.  2,  No.  1,  pp.  94-104,  1989. 

6)  Kleiser,  L„  and  T.  A.  Zang,  "Numerical  Simulation  of  Transition  in  Wall-Bounded  Shear  Flows", 
to  appear  Ann.  Rev.  Fluid  Mech.,  Vol.  23,  1991. 

7)  Mueller,  T.  J.,  "The  Role  of  Smoke  Visualization  and  Hot-Wire  Anemometry  in  the  Study  of 
Transition",  Recent  Advances  in  Turbulence,  edited  by  R.  E.  A.  Arndt  and  W.  K.  George,  Hemi¬ 
sphere  Press,  Washington,  1987. 

8)  Herbert,  T„  "Secondary  Instability  of  Boundary  Layers",  Ann.  Rev.  Fluid  Mech.,  Vol.  20,  pp. 
487-526,  1988. 

9)  Herbert,  T.,  F.  Bertolotti,  and  G.  Santos,  "Floquet  Analysis  of  Secondary  Instabilities  in  Shear 
Flows",  Instability  of  Time  Dependent  and  Spatially  Varying  Flows,  edited  by  D.  L.  Dwoyer  and 
M.  Y.  Hussaini,  Springer- Verlag,  pp.  43-57,  1985. 

10)  Kcgelman,  J.  T.,  and  T.  J.  Mueller,  "Experimental  Studies  of  Spontaneous  and  Forced  Transition 
on  an  Axisymmctric  Body",  AJAA  Paper  84-0008,  presented  at  AIAA  22nd  Aerospace  Sciences 
Meeting,  Reno,  Nevada,  January  9-12,  1984. 

11)  Kao,  K.-H.,  and  C.-Y.  Chow,  "Stability  Analyses  of  Boundary  Layer  on  a  Semi-Infinite  Circular 
Cylinder  With  and  Without  Spin",  AJAA  Paper  90-0116,  presented  at  28th  Aerospace  Sciences 
Meeting,  Reno,  Nevada,  January  8-11,  1990  . 

12)  Seiff,  A.,  and  B.  J.  Short,  "An  Investigation  of  Supersonic  Turbulent  Boundary  Layers  on  Slender 
Bodies  of  Revolution  in  Free  Flight  by  Use  of  a  Mach-Zehnder  Interferometer  and  Shadow¬ 
graphs",  NACA  TN-4364,  1958. 

13)  Whitfield,  J.  D.,  and  F.  A.  Iannuzzi,  "Experiments  on  Roughness  Effects  on  Cone  Boundary-Layer 
Transition  Up  to  Mach  16",  AJAA  J.,  Vol.  7,  No.  3,  pp.  465-470,  1969. 

14)  Stetson,  K.  F.,  Thompson,  E.  R.,  Donaldson,  J.  C.,  and  Siler,  L.  G.,  "Laminar  Boundary  Layer 
Stability  Experiments  on  a  Cone  at  Mach  8,  Part  1:  Sharp  Cone",  AJAA  Paper  83-1761,  July 
1983. 

15)  Fischer,  M.  C.,  and  L.  M.  Weinstein,  "Cone  Transitional  Boundary-Layer  Structure  at  M,-\4", 
AJAA  J.,  Vol.  10,  No.  5,  pp.  669-701,  1972. 

16)  Ericsson,  L.  E.,  "Effect  of  Nose  Bluntness  and  Cone  Angle  on  SIcnder-Vehicle  Transition",  AJAA 
J.,  Vol.  26,  No.  10,  pp.  1168-1174,  1988. 

17)  Pate,  S.  R.,  "Measurements  and  Correlations  of  Transition  Reynolds  Numbers  on  Sharp  Slender 
Cones  at  High  Speeds",  AJAA  J.,  Vol.  9,  No.  6,  pp.  1082-1089,  1971. 

18)  Wright,  R.  L.,  and  E.  V,  Zoby,  "Flight  Boundary  Layer  Transition  Measurements  on  a  Slender 
Cone  at  Mach  20",  AIAA  Paper  77-719  presented  at  10th  AIAA  Fluid  and  Plasma  Dynamics 
Conference,  Albuquerque,  New  Mexico,  June  27-29,  1977. 

19)  Mack,  L.  M.,  "Stability  of  Axisymmctric  Boundary  Layers  on  Sharp  Cones  at  Hypersonic  Mach 
Numbers",  AJAA  Paper  87-1413,  presented  at  AIAA  19th  Fluid  Dynamics,  Plasma  Dynamics  and 
Lasers  Conference,  Honolulu,  Hawaii,  June  8-10,  1987. 


-  28  - 


20)  Duck,  P.  W.,  "The  Inviscid  Axisymmctric  Stability  of  the  Supersonic  Flow  Along  a  Cylinder",  J. 
Fluid  Mech.,  Vol.  214,  pp.  61 1-637,  1990. 

21)  Malik,  M.  R.,  and  R.  E.  Spall,  "On  the  Stability  of  Compressible  Flow  Past  Axisymmetric 
Bodies",  High  Technology  Report  No.  1ITC-8905,  Hampton,  Virginia,  October  1989. 

22)  Macaraeg,  M.  G.,  and  Q.  I.  Daudpota,  "Numerical  Studies  of  Transverse  Curvature  Effects  on 
Compressible  Flow  Stability",  submitted  to  Phys.  Fluids. 

23)  Pruett,  C.  D.,  and  C.  L.  Street!,  "A  Spectral  Collocation  Method  for  Compressible,  Nonsimilar 
Boundary’  Layers",  submitted  to  lnt.  J.  Num.  Meth.  Fluids. 

24)  Ng,  L„  and  G.  Erlebacher,  "Secondary  Instabilities  in  Compressible  Boundary  Layers",  ICASE 
Report  No.  90-58  (NASA  Contractor  Report  182091),  1990. 

25)  Thompson,  P.  A.,  Compressible  Fluid  Dynamics,  McGraw-Hill,  1972,  pp.  632-635. 

26)  White,  F.  M„  Viscous  Fluid  Flow,  McGraw-Hill,  1974,  pp.  347-350. 

27)  Vachal,  J.,  "Phase  Ilia  Technical  Review-Environmental  Studies  and  Innovative  Concepts", 
NASA  High-Speed  Civil  Transport  Studies,  The  Boeing  Company,  Oct.  5,  1989. 

28)  "Phase  Ilia  High  Speed  Civil  Transport  Status  Review"  (NASA  Contract  NAS  1-18378),  McDon¬ 
nell  Douglas  Corporation,  Oct.  5,  1989. 

29)  Fetterman,  D.  E.,  Jr.,  "Preliminary  Sizing  and  Performance  of  Aircraft",  NASA  TM-86357,  July 
1985. 

30)  Gaster,  M.,  "A  Note  on  the  Relation  between  Temporally-Increasing  and  Spatially-Increasing  Dis¬ 
turbances  in  Hydrodynamic  Stability",/.  Fluid  Mech.,  Vol.  14,  pp.  222-224,  1962. 

31)  Mack,  L.  M.,  "Boundary-Layer  Linear  Stability  Theory",  Special  Course  on  Stability  and  Transi¬ 
tion  of  Laminar  Flow,  edited  by  R.  Michel,  AGARD  Report  No.  709,  pp.  3-1,3-81,  1984.  1990. 

32)  Macaraeg,  M.  G.,  C.  L.  Streett,  and  M.  Y.  Hussaini,  "A  Spectral  Collocation  Solution  to  the 
Compressible  Stability  Eigenvalue  Problem",  NASA  TP -2858,  1988. 

33)  Malik,  M.  R.,  "Prediction  and  Control  of  Transition  in  Supersonic  and  Hypersonic  Boundary 
Layers",  AIAA  J.,  Vol.  27,  No.  1 1,  1989. 

34)  Morkovin,  M.  V.,  "Transition  at  Hypersonic  Speeds",  ICASE  Interim  Report  1  (NASA  Contractor 
Report  178315),  May  1987. 

35)  Williamson,  J.  H.,  "Low-Storage  Runge-Kulta  Schemes”,  /.  Comp.  Phys.,  Vol.  35,  pp.  48-56, 
1980. 


streamwise  velocity,  u*/u 


2)  Streamwise  evolution  of  mean  streamwise  velocity  and  mean  temperature  in  the  boundary  layer  of 
a  cylinder. 


Growth  of  Displacement  Thickness  with  Mach  No. 


3)  Growth  with  increasing  Mach  number  of  relative  boundary-layer  displacement  thickness  8*/L* 
a  flat  plate. 


Normalized  Displacement  Thickness  vs.  C 


on 


4)  Variation  of  relative  displacement  thickness  8 */L*  with  curvature  C  at  Mach  4.S  (wind  tunnel 
conditions). 
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«\  Neutral-stabilitv  curves  for  first  and  second-mode  disturbances  at  Mach  4.5,  with  rays  of  constant 

dimensional  wavenumber.  _ a*; _ cC; _ “i  (adapted  from  L.  Mack  [31], 

with  permission.) 


7)  Influence  of  moderate  transverse  curvature  on  axisymmetric  first  and  second-mode  disturbances  at 
Mach  4.5. 


Growth  Rate, 


Temporal  Growth  Rate  vs.  Transverse  Curvature 


Temporal  Growth  Rate  vs.  Transverse  Curvature 


Transverse  Curvature,  C 


Influence  of  increasing  transverse  curvature  C  on  second-mode  disturbances  at  Mach  4.5,  for 
wavenumbers:  a)  a=0.20;  and  b)  a?=0.22 
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primary  wave  amplitude,  £, 


9)  Growth  rates  of  a  secondary  disturbances  of  fundamental  and  subharmonic  types  as  a  functions  of 
primary  disturbance  amplitude  at  Mach  1.6. 


Secondary  Instability 


10)  Influence  of  detuning  parameter  h  on  secondary  disturbance  growth  rate  at  Mach  4.5. 
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11)  Growth  rate  of  a  subharmonic  secondary  disturbance  as  a  function  of  primary  disturbance  ampli¬ 
tude  at  Mach  4.5. 
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Transverse  Curvature,  C 

12)  Influence  of  increasing  transverse  curvature  C  on  subharmonic  secondary  instability  at  Mach  4.5. 
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13)  Eigenfunction  components  of  a  subharmonic-type  secondary  disturbance  at  Mach  4.5. 


Subharmonic  Instability 


14)  Influence  of  increasing  Reynolds  number  on  a  subharmonic  secondary  instability  mechanism  at 
Mach  4.5. 
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Temporal  Growth  of  Selected  Harmonics 


15)  Temporal  growth  of  the  energy  content  of  selected  harmonics:  —  direct  simulation,  primary;  - 
-  -  direct  simulation,  secondary;  »  theory,  primary;  O  theory,  secondary;  A  theory  modified 
according  to  eqs.  (20);  □  growth  at  the  rale  Im((i)i-t  ^), 
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16)  Evolution  of  the  temperature  component  of  the  fundamental  (1,0)  mode. 


Temporal  Evolution  of  Mode  (1/2,1) 


17)  Evolution  of  the  temperature  component  of  the  (1/2,1)  harmonic  mode, 
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19)  Comparison  between  direct  numerical  simulation 
Reynolds  stress. 
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20)  Approximate  spatial  evolution  of  the  xXiI  Reynolds  stress  based  on  Gaster’s  transformation:  a) 
secondary  instability  theory;  b)  direct  numerical  simulation. 


wall- normal  coordinate 


Report  Documentation  Page 

x\<e  •Vrr^*y<3»or 


1.  Report  No.  2.  Government  Accession  No. 

NASA  CR- 1875 38 
ICASE  Report  No.  91-30 


4.  Tine  dno  Subtitle 


ON  THE  NONLINEAR  STABILITY  OF  A  HIGH-SPEED, 
AXISYMHETRIC  BOUNDARY  LAYER 


3.  Recipient's  Catalog  No. 


5.  Report  Date 

March  1991 


6.  Performing  Organization  Code 


7.  Author!  s) 

C.  David  Pruett 
Lian  L.  Ng 
Gordon  Erlebacher 


9.  Performing  Organization  Name  and  Address 

Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 
Hampton,  VA  23665-5225 


12.  Sponsoring  Agency  Name  and  Address 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  VA  23665-5225 


8.  Performing  Organization  Report  No. 

91-30 


10  Work  Unit  No 

505-90-52-01 


11.  Contract  or  Gram  No. 

NAS1-18605 


13.  Type  of  Report  and  Period  Covered 

Contractor  Report 


14  Sponsoring  Agency  Code 


15.  Supplementary  Notes 

Langley  Technical  Monitor: 
Michael  F.  Card 


To  be  submitted  to  Physics  of  Fluids. 


Final  Report 


J  ' fiU  The  stability  of  a  high-speed,  axisymmetric  boundary  layer  is  investi¬ 
gated  using  secondary  instability  theory  and  direct  numerical  simulation.  Parametric 
studies  based  on  temporal  secondary  instability  theory  identify  subharmonic  secondary 
instability  as  a  likely  path  to  transition  on  a  cylinder  at  Mach  4.5.  The  theoreti¬ 
cal  predictions  are  validated  by  direct  numerical  simulation  of  temporally-evolving 
primary  and  secondary  disturbances  in  an  axisymmetric  boundary-layer  flow.  At  small 
amplitudes  of  the  secondary  disturbance,  predicted  growth  rates  agree  to  several  sig¬ 
nificant  digits  with  values  obtained  from  the  spectrally-accurate  solution  of  the 
compressible  Navier-Stokes  equations.  Qualitative  agreement  persists  to  large  ampli¬ 
tudes  of  the  secondary  disturbance.  Moderate  transverse  curvature  is  shown  to  sig¬ 
nificantly  affect  the  growth  rate  of  axisymmetric  "second  mode"  disturbances,  the 
likely  candidates  of  primary  instability.  The  influence  of  curvature  on  secondary 
instability  is  largely  indirect  but  most  probably  significant,  through  modulation  of 
the  primary  disturbance  amplitude.  Subharmonic  secondary  instability  is  shown  to  be 
predominantly  inviscid  in  nature,  and  to  account  for  spikes  in  the  Reynolds  stress 
components  at  or  near  the  critical  layer. 


17.  Key  Word*  (Suggested  by  Author(sl) 

transition,  stability,  secondary  instabil¬ 
ity,  cylinder,  nonlinear,  direct  numeri¬ 
cal  simulation,  compressible,  Navier- 
Stokes 


19  Security  Clessif  (of  tfn*  report) 

Unclassified 


18.  Distribution  Statement 

02  -  Aerodynamics 

34  -  Fluid  Mechanics  and  Heat  Transfer 


Unclassified  -  Unlimited 


20  Security  Clissif  (of  this  psgsl 

2V  No.  of  pages 

Unclassified 

46 

22  Pries 

A0  3 


NASA-Langley,  1991 


